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Abstract 



New, gauge-independent, second-order Lagrangian for the motion of classical, 
charged test particles is used to derive the corresponding Hamiltonian formulation. 
For this purpose a Hamiltonian description of theories derived from the second-order 
Lagrangian is presented. Unlike in the standard approach, the canonical momenta 
arising here are explicitely gauge-invariant and have a clear physical interpretation. 
The reduced symplectic form obtained this way is equivalent to Souriau's form. 
This approach illustrates a new method of deriving equations of motion from field 
equation. 

1 Introduction 

In [Q a new method of deriving equations of motion from field equations was proposed. 
The method is based on an analysis of the geometric structure of generators of the Poincare 
group and may by applied to any special-relativistic, lagrangian field theory. In the case 
of classical electrodynamics, this method leads uniquely to a manifestly gauge- invariant, 
second order Lagrangian £ for the motion of charged test particles: 



C = L partide + C int = -Vl-v 2 (m - aVM^(t, Q, v 
where w M denotes the (normalized) four-velocity vector 



)) 



(1) 




(2) 
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and := u u V u u ,J ' is the particle's acceleration (we use the Heaviside-Lorentz system 
of units with the velocity of light c = 1). The skew-symmetric tensor MJ£*(t, q, v) is 
equal to the amount of the angular-momentum of the field, which is acquired by our 
physical system, when the Coulomb field accompanying the particle moving with velocity 
v through the space-time point (£, q), is added to the background (external) field. More 
precisely: the total energy-momentum tensor corresponding to the sum of the background 
field j^y and the above Coulomb field decomposes in a natural way into a sum of 1) terms 
quadratic in the background field, 2) terms quadratic in the Coulomb field 3) mixed terms. 
The quantity M^™* is equal to this part of the total angular-momentum M M „, which we 
obtain integrating only the mixed terms of the energy-momentum tensor. 

The above result is a by-product of a consistent theory of interacting particles and 
fields (cf. [||, ||), called Electrodynamics of Moving Particles. 

We have proved in |]J that the new Lagrangian (|I|) differs from the standard one 

L = L partide + L int = -Vl - v 2 (m - eu^A^t, q)) , (3) 

by (gauge-dependent) boundary corrections only. Therefore, both Lagrangians generate 
the same equations of motion for test particles in an external field. In the present pa- 
per we explicitly derive these equations and construct the gauge-invariant Hamiltonian 
description of this theory. 

Standard Hamiltonian formalism, based on the gauge-dependent Lagrangian (H), leads 
to the gauge-dependent Hamiltonian 



Hit, q, p) = y/m 2 + (p + eA(t, q)) 2 + eA (t, q) , (4) 

where the gauge- dependent quantity 

Pk ■= Pk n ~ eA k {t, q) = mu k - eA k (t, q) (5) 

plays role of the momentum canonically conjugate to the particle's position q k . 

As was observed by Souriau (see ||), we may replace the above non-physical momen- 
tum in the description of the phase space of this theory by the gauge-invariant quantity 
p kin . The price we pay for this change is, that the canonical contact form, corresponding 
to the theory of free particles: 

Q = dpf 1 A dq» , (6) 

has to be replaced by its deformation: 

n s :=n-ef^dq»Adq\ (7) 

where e is the particle's charge. 

Both Q and Qs are defined on the "mass-shell" of the kinetic momentum, i. e. on 
the surface (p fcm ) 2 = — m 2 in the cotangent bundle T*M over the space-time M (we 
use the Minkowskian metric with the signature (—,+,+,+)). The forms contain the 
entire information about dynamics: for free particles the admissible trajectories are those, 
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whose tangent vectors belong to the degeneracy distribution of f2. Souriau noticed that 
replacing (|j) by its deformation (0) we obtain the theory of motion of the particle in a 
given electromagnetic field f^ v . 

The new approach, proposed in the present paper is based on Lagrangian ([[J). It leads 
directly to a perfectly gauge-invariant Hamiltonian, having a clear physical interpretation 
as the sum of two terms: 1) kinetic energy mw and 2) "interaction energy" equal to the 
ammount of field energy acquired by our physical system, when the particle's Coulomb 
field is added to the background field. 

When formulated in terms of contact geometry, our approach leads uniquely to a new 
form Qjy: 

Q N := Q-eh^dq** Adq u , (8) 

where 

V := 2 (iW ~ u^f u]x u y ) (9) 
(brackets denote antisymmetrization), i. e. we prove the following 

Theorem 1 The one dimensional degeneracies of the form Qn restricted to the particle's 
"mass-shell" correspond to the trajectories of a test particle moving in external electro- 
magnetic field. 

It is easy to see that both Q s and fl^, although different, have the same degeneracy 
vectors, because h and / give the same value on the velocity vector u v : 

u u h^ v = u v f^ . (10) 

Hence, both define the same equations of motion. We stress, however, that our l^jv is 
uniquely obtained from the gauge-invariant Lagrangian ([]]) via the Legendre transforma- 
tion. 

The paper is organized as follows. In section |] we sketch briefly the (relatively little 
known) Hamiltonian formulation of theories arising from the second order Lagrangian. 
In section ||| we prove explicitly that the Euler-Lagrange equations derived from C are 
equivalent to the Lorentz equations of motion. Finally, Section |] contains the gauge- 
invariant Hamiltonian structure of the theory. 



2 Canonical formalism for a 2-nd order Lagrangian 
theory 

Consider a theory described by the 2-nd order lagrangian L = L(q l , q 1 , if) (to simplify 
the notation we will skip the index "i" corresponding to different degrees of freedom q l ; 
extension of this approach to higher order Lagrangians is straightforward). Introducing 
auxiliary variables v = q we can treat our theory as a 1-st order one with lagrangian 
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constraints := q — v = on the space of lagrangian variables (q,q,v,v). Dynamics is 
generated by the following relation: 

dL(q,v,v) = -^(pdq + n dv) = pdq + pdq + it dv + it dv . (11) 

where (p, it) are momenta canonically conjugate to q and v respectively. Because L is 
defined only on the constraint submanifold, its derivative dL is not uniquely defined and 
has to be understood as a collection of all the covectors which are compatible with the 
derivative of the function along constraints. This means that the left hand side is defined 
up to p,(q — v), where /i are Lagrange multipliers corresponding to constraints = 0. 
We conclude that p = \x is an arbitrary covector and ([TT]) is equivalent to the system of 
dynamical equations: 

dL 

di) 
dL 

^ dq 

dL , . 

»=&-'■ (12) 

The last equation implies the definition of the canonical momentum p: 

_dL . _dL d_ (dL\ 
^ dv dv dt \ dv J 

We conclude, that equation 

. d (dL\ d 2 (dL\ 



(15) 



dt 2 

is equivalent, indeed, to the Euler-Lagrange equation: 
5L d 2 (dL\ d (dL\ dL „ 





(-) 


-*\ 


(-) 


dt 2 


[dv / 


dt 


[dv J 



5q dt 2 \dv I dt \dv I dq 



The hamiltonian description (see e. g. is obtained from the Legendre transformation 
applied to (|TT|): 



— dH = pdq — qdp + 7r dv — v dn , (16) 

where H(q,p, v,it) = pv + n v — L(q, v, v). In this formula we have to insert v = v(q, v, tt), 
calculated from equation it — jk. Let us observe that H is linear with respect to the 
momentum p. This is a characteristic feature of the 2-nd order theory. 

In generic situation, Euler-Lagrange equations ( |I5] ) are of 4-th order. The correspond- 
ing 4 hamiltonian equations describe, therefore, the evolution of q and its derivatives up 
to third order. Due to Hamiltonian equations implied by relation (0), the information 
about succesive derivatives of q is carried by (v,ir,p): 
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v describes q 

dH , . 

4=-^. (17) 

hence, the constraint (p = is reproduced due to linearity of H with respect to p, 
tt contains information about q: 

dH , . 

*=ftF- (18) 

p contains information about 9 
dH dL 
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ir-P> ( 19 ) 



dv dv 

• the true dynamical equation equals 

. _dH_ _ dL 

^ dq dq 

3 Equations of motion from the variational principle 

In this section we explicitly derive the particle's equations of motion from the variational 
principle based on the gauge- invariant Lagrangian ([!]). The Euler-Lagrange equations for 
a second order Lagrangian theory are given by 

where, as we have seen in the previous section, the momentum pk canonically conjugate 
to the particle's position q k is defined as: 

dC . , 

Pk '- = d^~ * k ( } 

and 

dC 1 

^ : =^ = ^T^^ M r^,q,v). (23) 

Now, 

u"M%* = u°M% + u l Mtf = -u°r£* + u l e M m ^ , (24) 

where r^ n * and s™* are the static momentum and the angular momentum of the interaction 
tensor. They are defined as follows: we consider the sum of the (given) background field 
ffat and the boosted Coulomb field accompanying the particle moving with constant 
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four- velocity u and passing through the space-time point y = (t, q). Being bi-linear in 
fields, the energy-momentum tensor T total of the total field 

c tai ■= u» + (25) 

may be decomposed into three terms: the energy-momentum tensor of the background 
field T^ ield , the Coulomb energy-momentum tensor T partlcle , which is composed of terms 
quadratic in ffy u ^ and the "interaction tensor" T mf , containing mixed terms: 

rptotal rpfield _|_ rpparticle _|_ rpint (26) 

Interaction quantities (labelled with "int" ) are those obtained by integrating appropriate 
components of T mt . Because all the three tensors are conserved outside of the sources 
(i. e. outside of two trajectories: the actual trajectory of our particle and the straight 
line passing through the space-time point y with four- velocity u), the integration gives 
the same result when performed over any asymptoticaly flat Cauchy 3-surface passing 
through y. 

In particular, r mt and s int may be written in terms of the laboratory-frame components 
of the electric and magnetic fields as follows: 

rf(t,q,v) = / d 3 x (x k -q k )(BB + BB ) , (27) 

s%% q, v) = e mij ^ d 3 x (x* - <f )(D x B + D x B) J , (28) 

where D and B are components of the external field /, whereas D and B are components 
of f i.e.: 



v 2 



D ° (X;q ' V) = 4.|x-q|3 A , 2 , v(x _ q) ^ 3/2(X " q) ' ^ 



1-V 2 , , 

|x-q| 

B (x;q,v) = vxD (x;q,v). (30) 

It may be easily seen that quantities r k nt and s 1 ^ are not independent. They fulfill 
the following condition: 

= . (31) 

To prove this relation let us observe that in the particle's rest-frame (see the Appendix 
for the definition) the angular momentum corresponding to T mt vanishes (cf. |l|). When 
translated to the language of laboratory frame, this is precisely equivalent to the above 
relation. 

Inserting ([JTJ) into (|24]) we finally get 

*> = - I s '- + r^) • ' 32 > 
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The quantity r\ nt depends upon time via the time dependence of the external fields 
(D(t, x), B(t, x)), the particle's position q and the particle's velocity v, contained in 
formulae (|29|) - (0) for the particle's Coulomb field. 
Now, we are ready to compute p k from (|2^): 

mvk ^diii ( dTi k tdiTk , dir k \ 

~ Pk \dt +V dql ) V [drf dv>) ■ [ ' 

Observe, that the momentum p k depends upon time, particle's position and velocity but 
also on particle's acceleration. Hovewer, using (|27D one easily shows that due to the 

Lemma 1 

d7Xk 9711 0, (34) 



dv l dv k 

the term proportional to v vanishes (see Appendix for the proof). Moreover, one can 
prove the following 



Lemma 2 



dt dq l 



where we denote 

pl nt (t, q, v) = / d 3 x (D x B + D x B) k . (36) 

For the proof see Appendix. We see that p l k nt is the spatial part of the "interaction 
momentum" : 

pf(f,q,v) = Jj^dZ" , (37) 

where S is any hypersurface intersecting the particle's trajectory at the point (t,q(t)). 
The above integral is well defined (cf. 0) and it is invariant with respect to changes of 
E, provided the intersection point with the trajectory does not change. It was shown in 
[|l] that p™* is orthogonal to the particle's four- velocity, i.e. p l ^ l u^ = 0. 

Finally, the momentum canonically conjugate to the particle's position equals: 

P*=rf in +P* rf (t,q,v) • (38) 

It is a sum of two terms: kinetic momentum p^ n and the amount of momentum p l k nt which 
is acquired by our system, when the particle's Coulomb field is added to the background 
(external) field. We stress, that contrary to the standard formulation based on (|3|), our 
canonical momentum (1381) is gauge- invariant. 
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Now, Euler-Lagrange equations fl2~T| ) read 

^ + ^ = ^ (39) 



dpf n dp{ nt dC 
dt dt dq k 



or in a more transparent way: 

d ( mv k \ f dpf* jdpf*} rf dpf* dirA 

dt\^T^J { 8t + dq l J \dv l dq k ) ' { } 

Again, using definitions of iii and p k nt one shows that due to the following 



Lemma 3 



. (41) 



dv l dq k 

the term proportional to the particle's acceleration vanishes (for the proof see Appendix). 

The last step in our derivation is to calculate — h i>'-^r- In the Appendix we show 
that the following identities hold: 



Lemma 4 



^ + vl ^f = ~ e uUfk ^ q) = ~< Ek ^ «l) + £kimV l B m (t, q)) . (42) 



Therefore, the term — \-v q~t gives exactly the Lorentz force acting on a test particle. 

This way we proved that the Euler-Lagrange equations (|21~D for the variational problem 
based on £ are equivalent to the Lorentz equations for the motion of charged particles: 



d ( mt)(. 



dt \ 



e(E k (t, q) + e klm v l B m (t, q)) . (43) 



4 Hamiltonian formulation 



By Hamiltonian formulation of the theory we understand, usually, the phase space of 
Hamiltonian variables V = (q,p) endowed with the symplectic 2-form u = dp A dq and 
the Hamilton function H (Hamiltonian) defined on V . This function is interpreted as 
an energy of the system. However, for time-dependent systems this framework is usually 
replaced by (a slightly more natural) formulation in terms of a contact form. For this 
purpose one considers the evolution space PxR endowed with the contact 2-form (i.e. 
closed 2-form of maximal rank): 

u H := dp Adq- dH Adt . (44) 

In analytical mechanics this form, or rather its "potential" pdq — Hdt, is called the 
Poincare-Cartan invariant. Obviously, ujh is degenerate onPxRand the one- dimensional 
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characteristic bundle of ujh consists of the integral curves of the system inPxR. This 
kind of description may be called the "Heisenberg picture" of classical mechanics. In this 
picture states are not points in V but "particle's histories" inPxR (see [|J). 

Let us construct the Hamiltonian structure for the theory based on our second order 
Lagrangian C Let V denote the space of Hamiltonian variables, i.e. (q, p, v,7r), where 
p and 7r stand for the momenta canonically conjugate to q and v respectively. Since our 
system is manifestly time-dependent (via the time dependence of the external field) we 
pass to the evolution space endowed with the contact 2-form 

■= dp k A dq k + dn k A dv k — dH A dt , (45) 

where 7i denotes the time-dependent particle's Hamiltonian. 

To find Ti on V x R one has to perform the (time-dependent) Legendre transformation 
(q, q, v,v) — > (q, p, v,7r), i.e. one has to calculate q and v in terms of Hamiltonian 
variables from formulae: 

p*=a^-**> ** = w (46) 

This transformation is singular due to linear dependence of £ on v and gives rise to the 
time-dependent constraints, given by equations ( |2"BD and (^). The constraints can be 
easily solved i.e. momenta p k and TT k can be uniquely parameterized by the particle's 
position q k , velocity v k and the time t. Let V* denote the constrained submanifold of 
the evolution space PxR parametrized by (q, p km ,t). The reduced Hamiltonian on V* 
reads: 

777 

H(t, q, v) = p k v k + n k v k -£= + v k p k nt (t, q, v). (47) 

V 1 — 

Due to identity u^p™ 1 = (cf. [|IJ) we have 

777 

H(t, q, v) = -== - p™% q, v) , (48) 



v 



and, therefore, 



Theorem 2 The particle 's Hamiltonian equals to the u —po " component of the following, 
perfectly gauge-invariant, four-vector 

P, ■= P h ; n + Pt% q, v) = mu, + pf% q, v). (49) 
Using the laboratory-frame components of the external electromagnetic field we get: 

p**(t, q, v) = - J d 3 x (DD + BB ) . (50) 

Now, let us reduce the contact 2-form (f45|) on V* . Calculating p k = p k ((i,p km ,t) and 
TTfe = ^(q? p km , t) from (23) - (|38|) and inserting them into (^5|) one obtains after a simple 
algebra: 

n N = dp k J n Adq^-e dq» A dq u , (51) 
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where q° = t and is the following 4- dimensional tensor: 

faint faint 

eM<,q,v):=^-^. (52) 
Using techniques presented in the Appendix one easily proves 
Lemma 5 



dp" 



int 



erL A A, , (53) 



where 



n M A := (5/ + «X (54) 

is the projection on the hyperplane orthogonal to w M (i.e. to the particle's rest-frame 
hyperplane, see the Appendix). Therefore 

V = ii m a /a, - n„ A / AAt = 2(/^ - u [fl f u]x u x ) . (55) 

where a^ft^j := \{a a bp — a^b a ). The form Q N is defined on a submanifold of cotangent 
bundle T*M defined by the particle's "mass shell" (p fem ) 2 = —m 2 . 

Observe, that the 2-form (|5l]) has the same structure as the Souriau's 2- form ([?]). 
They differ by the "curvature" 2-forms / and h only. However, the difference u h — /" 
vanishes identically along the particle's trajectories due to the fact that both f^ v and 
have the same projections in the direction of (see formula fllCf)). We conclude that the 
characteristic bundle of Qn and Qs are the same and they are described by the following 
equations: 

q k = v k , (56) 



v 



VT^ -(g kl - v k v l ){E t + eujv'Bi) , (57) 



7/i 



which are equivalent to the Lorentz equations 

We have two different contact structures which have the same characteristic bundles. 
Therefore, from the physical point of view, these forms are completely equivalent. 



Appendix 

Due to the complicated dependence of the Coulomb field D and B on the particle's 
position q and velocity v, formulae containing the respective derivatives of these fields 
are rather complex. To simplify the proofs, we shall use for calculations the particle's 
rest-frame, instead of the laboratory frame. The frame associated with a particle moving 
along a trajectory ( may be defined as follows (cf. [|J, at each point (t, q(t)) G C we 
take the 3-dimensional hyperplane S 4 orthogonal to the four-velocity (the rest-frame 
hyper surface) . We parametrize S 4 by cartesian coordinates (x k ), k = 1,2,3, centered at 
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the particle's position (i.e. the point x k = belongs always to (). Obviously, there are 
infinitely many such coordinate systems on Tit, which differ from each other by an 0(3)- 
rotation. To fix uniquely coordinates (x k ), we choose the unique boost transformation 
relating the laboratory time axis d/dy° with the four- velocity vector U := u**-^. Next, 
we define the position of the d/dx k - axis on E t by transforming the corresponding d/dy k 
- axis of the laboratory frame by the same boost. The final formula relating Minkowskian 
coordinates with the new parameters (t,x k ) may be easily calculated (see e. g. f3[) 
from the above definition: 

y°(t,x l ) := t + = x l Vi(t) , 

y/l ~ v 2 (0 

y k (t,x l ) := q k (t) + (6 k + ^ 2 )v k v l )x l , (A.l) 

l 



where we denote <p(z) :— , — ... , -. — , — . 

' v ' z Vvl— z I vl- z(l+vl— z) 

Observe, that the particle's Coulomb field has in this co-moving frame extremely 
simple form: 

^o(x) = i - J , B (x)=0, (A.2) 

where r := |x|. That is why the calculations in this frame are much easier than in the 
laboratory one. 

Let T>k and Bk denote the rest-frame components of the electric and magnetic field. 
They are related to Dk and Bk as follows: 

1 



v 



P fc (x,t;q,v) = _ _ K - VT^p(v 2 )v l v k ) D,(y) - e^B^y) , (A.3) 



B fc (x,t;q,v) = -, = (5 l k - VY^^v 2 )v l v k ) B t (y) + e^D^y) , (A.4) 



(the matrix (5 l k — \/l — v 2 Lp(v 2 )v l Vk) comes from the boost transformation). 

The field evolution with respect to the above non inertial frame is a superposition of 
the following three transformations (cf. [T|], [0, ||): 

• time-translation in the direction of U, 

• boost in the direction of the particle's acceleration a k , 

• purely spatial 0(3)-rotation around the vector u m , 
where 

J: 



a 



(*? + <^VW) v l , (A.5) 



u m := v(v 2 )v k v l e k i m . (A.6) 

V 1 — v 2 
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Therefore, the Maxwell equations read (cf. 

^ d 



V n = vT 



B n = Vl^v 1 



dx m 
d 

dx m 



(e mk p n - e nk J) m ) UkX 1 - e mn k {l + a'a^B* 
(e m *£" - t n \B m ) ui k x l + e mn fc (l + a i x l )V k 



(A.7) 
(A.8) 



(the factor yl — v 2 is necessary, because the time t, which we used to parametrize the 
particle's trajectory, is not a proper time along ( but the laboratory time). 

On the other hand, the time derivative with respect to the co-moving frame may be 
written as 



d _ d k d . k d ( d\ . k d 



dv h 



Therefore, taking into account (|A.7 ) and (|A.8 ) we obtain: 



vdt 



B" 



e nmk d m B k , 
-VT^ e nmk d m V k 



and 



u 



dv 1 
d 

— B n 

dv 1 



VT^2 Q 



du 



vT^2 9 



dv 

duJk 
dv 1 



_ ( / m k T} n 



n nk<T}m\ i " r mn r>k 

1 ) ~dv J k 



(A.9) 

(A.10) 
(A.ll) 

(A.12) 
(A.13) 



To calculate the derivatives of D k and B k with respect to the particle's position observe, 
that 



d 



dy k y/1 - v 
Therefore 



_d_ 

dx i 



_d_ 

dq k 
d 
dq k 



B n 



v k 



e nmi d m B i +(S i k + l p(v 2 )v i v h )d i V r 



Vk 



VT 



e nmi d m Vi + [51 + ip(v')v l v k ) diW 1 . 



(A.14) 

(A.15) 
(A.16) 



Now, using ( IA.10I ) (|A.13| ) and ( |A.15|) -( |A~T^ ) we prove Lemmas |T]-g 

1. Proof of Lemma [l]: 

Observe, that "interaction static moment" in the particle's rest frame reads: 



-pint 
K k 



x k (D T> + B B) d 3 x 



Att js* r 3 



XkX V'd 3 x 



(A.17) 
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Taking into account that 

we obtain the formula for n k in terms of Rf 1 *: 
1 



TTfc 



Now, using ( |A.12| ) one gets: 
d 



(5i + v (v 2 y Vk ) rt 1 



da 



v 



dcu r 



dv 



l --im 



3 J^int 



where 



X,; 



e 

An' 



A _ ^ijk 



X J X 



^B k d 3 x. 



Therefore 

dv 1 dv k 



Ai -pint 
A kl n i 



where 



.4' 



ki 



_d_ 

dv k 



+ e ■ 

1 jm 



Dim 



VT 

9a 



Dim v 



CM 



2 Vw fe 



(<5/ + ip(v 2 )v^vi 



dto m 
dv k 



(si + (p(v 2 y Vk 

9s / da 1 da 



dv 



W - {si + ^(v 2 ) 



V Vi 



da m 
dv k 



1 



9a* 9a" 



Using the following properties of the function (p(z): 

2ip(z)-(l-z)- 1 + z<p 2 (z) = 0, 
2( f /(z)-(l-z)- 1 <p(z)-<f?(z) = 0, 



(A.18) 



(A.19) 



(A.20) 



(A.21) 



(A.22) 



(A.23) 



(A.24) 



(A.25) 
(A.26) 



one easily shows that A l kl = 0. Moreover, observe that B k Y defined in ( A.24 ) is antisym- 
metric in (im). Therefore, to prove (p4]) it is sufficient to show that the quantity Xj m is 
symmetric in (im). Taking into account that B k = e klm diA m , where A m stands for the 
rest-frame components of vector potential, one immediatelly gets: 



-ijk 



X J X 



m B k d 3 x 



T (A k X )(3XiX m T Cji m ) d X , 



(A.27) 
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which ends the proof of 
2. Proof of Lemma |2|: 



To prove (|35|) observe that 



diik i d-Kk 



dt 



dq l 



dt 



at 



R 



it it 



(A.28) 



Now, using ( |A.10| ) we obtain 



R 



it it 



vT 



at 



Due to the Gauss theorem 



— e k ^ m d j B m d 3 x 



I a, ( ^e kjm B m ) d 3 x + I e km ^B m d 3 x 



km X k 



2 Q 



St 



^e kjm B m ) d 3 x 



XjXiXfc f c j r 



9s t r 



, e^ m B m da = 0, 



(A.29) 



(A.30) 



where da denotes the surface measure on <9£+. Moreover, observe that "interaction mo- 



mentum" in the particle's rest-frame reads: 

k 

P} nt :=e ikm f (V k B m + V k B™)d 3 x = ^e lkm [ X -B m d 3 x 
Js t 47T Js t r 4 

Therefore 

( —) Rip-t = y/l _ v 2 pint 

\dt) u 1 v 1 

Using the relation between p l k nt and P- nt 
rf" = (<5 + ¥>(vV«*) Pt nt 



(A.31) 



(A.32) 



(A.33) 



we finally get (|35|). 

3. Proof of Lemma |3|: 

Using ( |A.13| ) and (]A.15 ) we obtain: 



a 

pint 

dv l 1 

a 

" T}int 



vT 



du™ Jpmt _da™ 



(9?) 



Vl 



pint + (gj + v (v 2 )v j Vl ) Y 



where 



Y 



.1 y 



47T js t r 5 



(A.34) 
(A.35) 

(A.36) 
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Now, taking into account ( |A.19| ) and ( |A.33j ) we have 



dp l k d%i 
dv l dq k 

where 



C\ t Pr , (A.37) 



c\i ■= w (4 + ^ 2 y Vk ) - vr^ (si + ^y Vk ) — e jn : 

+ ( 5 i + V^ 2 )v l vi) ■ (A.38) 



One easily shows that due to properties ( |A.25| ) ( p\.26[ ) C l kl = 0, which ends the proof of 

Hi. 



4. Proof of Lemma 

Finally, to prove (fE3) let us observe that 



¥ + *w - (I) B # - w + «w>) (I) „ p ™' • < A ' 39 > 



Now, due to ( | A. 11| ) we get 



dt) u 1 A-KJdT.tr 2 



= -VI - v 2 — lim / —Vi da = -Vl - v 2 eVJt, 0) , (A.40) 

ATT r ^0j S (r ) r 2 

where we choose as two pieces of a boundary <9£ t a sphere at infinity and a sphere S(ro). 
Using the fact that in the Heaviside-Lorentz system of units T> k = E\~ and taking into 
account the formula ( |A.3| ) we finally obtain 
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